function [rID,cID,mat]=sta(Cls)
%%
%rCls/cCls : numSample*2
%
    %numSample = size(rCls,1);
    rCls = Cls(:,1);
    cCls = Cls(:,2);
    rID = sort(unique(rCls));
    cID = sort(unique(cCls));
    mat = zeros(length(rID),length(cID));
    for r = 1:length(rID)
        for c = 1:length(cID)
            idx = find(rCls == rID(r));
            num = length(find(cCls(idx) == cID(c)));
            mat(r,c) = num;
        end
    end
end
